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Abstract 

We suggest a new model for the dynamics of a suspension bridge through a system of nonlinear 
nonlocal hyperbolic differential equations. The equations are of second and fourth order in space 
and describe the behavior of the main components of the bridge: the deck, the sustaining cables 
and the connecting hangers. We perform a careful energy balance and we derive the equations from 
a variational principle. We then prove existence and uniqueness for the resulting problem. 


1 Introduction 

The main span of a suspension bridge is a complex structure composed by several interacting compo¬ 
nents, see Figure Four towers sustain two cables that, in turn, sustain the hangers. At their lower 
endpoint the hangers are linked to the roadway and sustain it from above. The roadway (or deck) 
may be seen as a thin rectangular plate. The hangers are hooked to the cables and the roadway is 
hooked to the hangers; the weight of the deck deforms the cable and stretches the hangers that exert 
a restoring action on the roadway. 



Figure 1: Sketch of a suspension bridge. 


This complex structure is extremely interesting from a mathematical point of view. It appears quite 
challenging to find a model simple enough to be mathematically tractable and sufficiently accurate to 
display the main features of a real bridge. One of the main problems of suspension bridges is their 
instability, in particular they are prone to torsional oscillations. This is apparent in videos available 
on the web |20j as well as in many other events, see e.g. |10j for a survey. Partial explanations of the 
instability are based on aerodynamic effects such as vortex shedding, flutter, parametric resonance, 
aerodynamic forces, see [6]. But none of these theories explains how a longitudinal oscillation can 
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rapidly switch into a torsional one. Scanlan m p.209] writes that the switch is due to some fortuitous 
condition, a justification which does not seem very scientific. And, indeed, McKenna |14) raises many 
doubts about the mathematical models and motivations used in classical literature. In a recent paper 
[2j we suggested that the onset of the torsional instability could be of purely structural nature. The 
model adopted there was fairly simplified and did not describe accurately the nonlinear behavior of the 
cables+hangers restoring force. By using Poincare-type maps we were able to show that a longitudinal 
oscillation in an ideally isolated bridge may switch rapidly into a torsional oscillation when enough 
energy is present within the structure. 

Aiming to display the same phenomenon of transfer of energy between longitudinal and torsional 
oscillations in an actual bridge, in this paper we suggest a new and realistic model which takes into 
account all the main components of the bridge and its quantitative structural parameters. We perform 
an energy balance and we derive the corresponding Euler-Lagrange equations by variational principles. 
The dynamics of the bridge is modeled through a system of four hyperbolic equations, three of them 
of second order in space and one of fourth order. This system also contains nonlinear couplings (the 
restoring action of the hangers) and nonlocal terms (the elongation of the extensible cables). In fact, 
the full energy balance is rather complicated, therefore we approximate the equations with a first order 
expansion and we prove that the initial-boundary value problem for the hyperbolic system admits a 
unique solution. In a forthcoming paper [3] we will study higher order approximations. 

In Section we derive an accurate model and we proceed in two steps. First, in Sections 2T and 


2.2 we compute the Lagrangian of a system composed by a one-dimensional beam suspended to a 
cable through hangers, see Figure]^ In the classical literature [71 113L [T6] one neglects either the mass 
of the cable, in which case the load is distributed per horizontal unit and the cable takes the shape 
of a parabola, or the mass of the beam, in which case the load is distributed per unit length and the 
cable takes the shape of a catenary. Since none of these masses is negligible, we take here into account 
both their contributions. This gives rise to a somehow intermediate shape and the linearized problem 


becomes a Sturm-Liouville problem with a weight. Second, in Section 2.3 we extend the result to the 
full suspension bridge by modelling the roadway as a degenerate plate, that is, a beam representing 
the midline of the roadway with cross sections which are free to rotate around the beam; simplified 
equations representing this model were previously analyzed in mm- The midline corresponds to 
the barycenter of the cross sections. The cross sections are seen as rods that can rotate around their 
barycenter. The endpoints of the cross sections are the edges of the plate that are connected to the 
sustaining cables through the hangers. 

Then we derive the Euler-Lagrange equations by variational methods, that is, by finding critical 
points of the Lagrangian. The resulting equations form a system of semilinear hyperbolic equations, 


three of them being of second order and one of fourth order, see (38). All the coefficients are continuous 


time-independent functions and (38) is a nonlocal differential problem due to the presence of the 
integral term representing the cables elongations. A nonlinear term is responsible for the coupling 
between the equations. These peculiarities show that (38) is not a standard problem and therefore we 


need to prove existence and uniqueness of a solution. This is done in Section 


2 The physical model 

Throughout this paper we denote the partial derivatives of a function y = y{x, t) by 


y = 


dx ’ 


y = 


dt 


and similarly for higher order derivatives. 
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2.1 A beam sustained by a cable through hangers 

As a first step, we derive the Lagrangian of a suspended beam connected to a sustaining cable through 
hangers. Let w be the position of the beam and p — she the position of the cable, as in Figure The 
purpose of this section is to derive the equation of this cable-hangers-beam system. 



Figure 2; Cable-beam structure modeling a suspension bridge. 


We model the cable as a perfectly flexible string subject to vertical loads. If we assume that the 
string has no resistance to bending, the only internal force is the tension F = F{x) of the string which 
acts tangentially to the position of the curve representing the string. In Figure [^we sketch a picture of 


Fcos6 



the string whose endpoints are A and B and whose position is described by a function —s(x) < 0, the 
origin being below B. The horizontal direction represents the abscissa x whereas the downwards axis 
represents the vertical displacement. Denote hy 9 = 9{x) the angle between the horizontal x-direction 
and the tangent to the curve so that 

— s'(x) = tan0(x). (1) 

The horizontal component of the tension is constant, that is, 

F{x)cos9{x) = Hq > h . (2) 

The classical deflection theory of suspension bridges models the bridge structure as a combination 
of a string (the sustaining cable) and a beam (the deck). The cable carries its own weight, the weight 
of the hangers (which we assume to be negligible) and the weight of the beam. Possible deformations 
of the cable and the beam are assumed to be small, that is, we aim to model 

small displacements of the beam and the cable. (3) 
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In Figure]^ we represent the position —s{x) of the cable at rest with —sq < 0 being the level of the 
left and right endpoints of the cable (so is the height of the towers). We denote by L the distance 
between the towers. Assume that a beam of length L and linear density of mass M is hanged to the 
cable whose linear density of mass is m: the segment OQ represents the horizontal position of the 
beam at rest. Since the spacing between hangers is small relative to the span, the hangers can be 
considered as a continuous membrane connecting the cable and the beam. Then the cable is subject 
to a downwards vertical force given by 

q{x) = + s'(x)^^ g 

g being the gravitational constant. In this situation, the vertical component of the tension has a 
variation given by 

— [F{x)sm6{x)] = —q{x) 

where we recall that the positive vertical axis is oriented downwards. In view of Q we then obtain 

Ho^[ta.ne{x)] = -q{x) 

Moreover, using 0 and taking into account that the cable is hanged to two towers of same height so, 
we get 



If one neglects the mass of the cable {m = 0) then the solution is the parabola 


Sp{x) = So - - x ), 

ZHO 

while if one neglects the mass of the beam (M = 0) the solution is the catenary 


(5) 


Sc{x) 


ih 

mg 


cosh 



(2x 



cosh 



+ So . 


We point out that if we take m = 0 in the above explicit equations and we assume that the sag-span 
ratio for is 1/12 (a typical value for actual bridges, see e.g. [IHl §15.17]), by using ([^ we have 

gM _ 2 
2J{^ ~ ^ 

and then 

0 < Sp{x) — Sc{x) < Sp{L/2) — Sc{L/2) « 6 • 10“^ L Vx G (0, L). 

Hence, for a typical length of the bridge L = 1km the maximum difference between the two configu¬ 
rations is about 6m. 

Our purpose is to consider both the masses of the beam and the cable (m, M > 0). Then Q does 
not have simple explicit solutions. Moreover, this second order equation has no variational structure 
and the problem has boundary conditions, therefore existence and uniqueness of the solution are not 
granted, but are proved in the next statement. 

Proposition 1. For any m, M > 0 there exists a unique solution s = s(x) of Q; this solution is 
symmetric with respect to x = L/2, that is, s{x) = s{L — x). 
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Proof. We first prove the existence and uniqueness of a symmetric solution of Q, then we show that 
any solution of is symmetric. These two steps will prove the statement. 

For any si G M the initial value problem 


Hqs”{x) = (^M + m^l + s'(x)2^ g , 

s{L/2) = Si , 

s'(L/2) = 0, 


( 6 ) 


has a unique solution defined on the whole real line for all values of si. Moreover, solutions of (|^ with 
different values of si differ by a constant and they are all symmetric with respect to x = L/2. Then 
there exists a unique si < sq such that the solution of Q satisfies s(0) = s{L) = sq: this solution is 
the unique symmetric solution of @- 

Take now a solution s of Q which we know to exist in view of what we just proved. Since s 
is convex and s(0) = s{L), by the Lagrange Theorem there exists a unique xq G (0,L) such that 
s'{xo) = 0. Hence, s solves the initial value problem 


Hqs''{x) = (^M + m^l + s'(x)2^ g , 
s(xo) = s(xo), 
s'(xo) = 0. 


(7) 


But also s(2xo — x) solves Q: then by uniqueness of the solution we infer that s{x) = s(2xo — x). 
Since s(0) = s{L) = sq and since s is convex, we finally deduce that 2xo = L. Therefore, any solution 
s of Q satisfies s(x) = s{L — x). □ 

Hence, for general m,M > 0 the cable takes the shape of a curve with a U-shaped graph, something 
in between a parabola and a catenary. The length at rest of the cable is 


Lc= -y/l + s'(x)^ dx . 


( 8 ) 


Let 


^(x) := Vl + s'(x)2 (9) 

be the local length of the cable, which is related to the unique solution of Q (see Proposition and 
will often appear in the computations. In classical literature this term is usually approximated with 
^(x) = 1. While modeling suspension bridges, Biot and von Karman [TJ p.277] warn the reader by 
writing ...whereas the defleetion of the beam may be considered small, the deflection of the string, i.e., 
the deviation of its shape from a straight line, has to be considered as of finite magnitude. Whence, Q 
is acceptable whereas no approximation of s'{x) should be taken. Nevertheless, at a subsequent stage, 
Biot-von Karman [3 (5.14)] decide to neglect s'(x)^ in comparison with unity, this leads to ^(x) = 1. 
As shown in |12j this mistake may significantly change the responses of the system. Needless to say, 
this rough approximation considerably simplifies all the computations. 


2.2 Lagrangian of the cable-hangers-beam system 

In this section we write all the components of the Lagrangian of the single cable-hangers-beam sys¬ 
tem. Since the energies involved are complicated, we will use Q in order to approximate them with 
asymptotic expansions up to second order terms. This will give rise to Euler-Lagrange equations with 
linear terms. 
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• Kinetic energy of the beam. If tc = w{x, t) denotes the position of the beam, then the kinetic 
energy of the beam is given by 


^ M r 2 , 

Ekb = ^ W dx. 
^ JO 


( 10 ) 


• Kinetic energy of the cable. If p{x,t) — s{x) denotes the position of the cable, then the kinetic 
energy of the cable is given by 


m 


Ekc= ^ / p ^{x)dx. 

^ Jo 


( 11 ) 


• Gravitational energy. If we neglect the small extensions of the cable, the gravitational energy is 
given by 


Ea = - 


g / (^Mw + mp^{x'^ dx . 


The minus sign is due to the downward orientation of the vertical axis. 

• Elastic energy in the hangers. The hangers behave as stiff linear springs when in tension and 
they do not apply any force when compressed. The latter case is called slackening of the hangers. 
Let A = A(x) be the length of the unloaded hanger at position x G (0, L); this is the vertical length 
of the hanger when it is fixed to the sustaining cable and before hanging the beam. After the beam 
is installed, the hangers are in tension and reach a new length s(a:) > A(x), where s is the solution 
of Q; if no additional load acts on the system, the equilibrium position of the beam is re (re) = 0 
(corresponding to the segment OQ in Figure and the position of the cable at equilibrium (beam 
installed) is —s(x). As one expects from a linear spring, the Hooke constant k(x) of each hanger is 
proportional to the inverse of its unloaded length, so that 


k(x) = 


kq 

A(x) 


(ko > 0). 


Then the (linear density of) force due to the hanger in position x is 

k(x) ^s(x) — A(x)^ . 

Since the (linear density of) weight of the beam is Mg, we find 


( 12 ) 


, . Ko(s(x) - A(x)) 

Mg = k(x)(^s(x) - A(x)J =-- 


(13) 


which shows that the relative elongation of the hangers is proportional to the density of weight after 
the beam is installed: this relationship readily gives the following dependence of the elongation s(x) 
on the unloaded length 

The actual length of the hanger is w(x,t) — (p(x,t) — s(x)), so that the hanger applies a force 

/ \ ^ 9 \ 

F{w — p) = k{x) yw{x, t) — {p{x, t) — s{x)) — A(x)j = k(x) yw{x, t) — p{x, t) H— j—j 


= (j^Ko + Mg) 


w{x, t) — p{x, t) 
s{x) 


+ 


MgY', 


(14) 


note that F{0) = Mg, that is, the restoring force of the hangers balances gravity at equilibrium. The 
positive part in (14) models the fact that the hangers behave as linear springs when extended and do 
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not yield any force when slackened. The equalities in (14) are due to (13). The elastic energy is then 
given by 


Eh{w-p) = ^ 


[w{x,t) -p{x,t) + 


dx 


k{x) (^w{x,t) — p{x,t) + s{x) — X{x)J j dx 


(15) 


Remark 1. The special form of F in (14) was first suggested by McKenna-Walter m and well 
models the possible slackening of the hangers, at least in a first order approximation as in the present 
paper. Brownjohn [9j claims that the slackening mechanism is not as simple as an on/off force, as 
described by the positive part of the elongation. Our results remain true if the force F is replaced by 
any Lipschitz continuous function. □ 

• Bending energy of the beam. The bending energy of a beam depends on its curvature, see [7] 
and also m for a more recent approach and further references. The Young modulus F is a measure of 
the stiffness of an elastic material and is defined as the ratio stress/strain. If Idx denotes the moment 
of inertia of a cross section of length dx, then the constant quantity FI is the flexural rigidity of the 
beam. The energy necessary to bend the beam is the square of the curvature times half the flexural 
rigidity: 

FI {w"Y 


Es=- 


; Y^l + {w'Y dx 


(16) 


/o (1 + (tc')^)^ 

where we highlighted the curvature and the arclength. This energy is not convex and it fails to be 
coercive in any reasonable functional space so that standard methods of calculus of variations do not 
apply. Assuming ([^, an asymptotic expansion in (16) yields 




'o + 2 , 


[w 


1 - 


/\2 




dx + o{w'^). (17) 


Therefore, from now on, we simply take 


eb = ^ [\wydx. 

^ Jo 


• Stretching energy of the beam. If large deformations are involved, the strain-displacement 
relation is not linear and a possible nonlinear model was suggested by Woinowsky-Krieger [23]: he 
modified the classical Bernoulli-Euler beam theory by assuming a nonlinear dependence of the axial 
strain on the deformation gradient and by taking into account the stretching of the beam due to its 
elongation in the longitudinal direction. In this situation there is a coupling between bending and 
stretching and the stretching energy is proportional to the elongation of the beam which results in 


Es = - 
2 


(yiTFF-i) 


dx 


7 


{y) dx 


where 7 > 0 is the elastic constant of the beam. Since this term is of fourth order, in the sequel we 
drop it. 

• Stretching energy of the cable. The tension of the cable consists of two parts, the tension at 
rest 

H{x) = Hoax) (18) 

and the additional tension 
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( 19 ) 


due to the increment of length r(p) of the cable. The latter requires the energy 

„ , , AE 

Etc{p) = ■ 

Assuming Q, we may approximate Etc{p) as follows. First, we note the asymptotic expansion 

Vl + [.'(X) - ^(x)Y - Vl + s'(xy = - + lly + o(p'(.^)^) . (20) 

Then we approximate the increment r(p) of the length of the cable, due to its displacement p, by 
using ( 1 ^: 


r(p) := / ^1 + [s'(x) -p'{x)Pdx - Lc 

Jo 


p'{x)‘^ s'{x)p'{x] 

2^(x)3 


e(x) 


By squaring we find 


Tipf = 


rL j 


s'{x)p'{x) 

i{x) 


\ 

dxj +o^(p')^^ 


Therefore, we approximate (19) by 


dx. 


= 2 l: 


rL j 


s'{x)p'{x) 

i{x) 


dx 


By taking the variation of this energy and integrating by parts we obtain 


dEtc{p)z = —— 
Lir 


rL J 


s'{x)p'{x) 

i{x) 


dx 


s"{x)z{x) ^ 

^ ^.3 dx + o{p) 


'0 


i{x) 


( 21 ) 


On the other hand, the amount of energy needed to deform the cable at rest under the tension 
(18) in the infinitesimal interval [x,x + dx] from the original position —s(x) to —s{x) + p{x) is the 
variation of length times the tension, that is 

E{x) dx = Ho^{x) ^-\/l + [s'(x) — p'(x)]2 — y^l + s'(x)2^ dx , 


hence, using ( 20 ) and neglecting the higher order terms, the energy necessary to deform the whole 
cable is 

p'{xY 


E{p) = j E{x)dx = Ho j s'(x)p'(x) + 


2^{xf 


dx. 


( 22 ) 


Then, the variation of energy a ^[E{p + az) — E{p)] as a —)• 0 and some integration by parts lead to 

'' p'{x)z'{x)' 


dE{p)z = Ho ^ - s'(x) 2 :'(x) + 


= Ho 


= Ho 




+ 


2s's''p' 


dx 


C(x)2 ^(x)4 

" 2s's''p' 




C(x)2 ^(x)4 


z{x) dx 
z(x) dx. 


(23) 


• The Euler-Lagrange equations. The Lagrangian of the system is obtained as the sum of the 
kinetic energies minus the sum of the potential energies. To compute the Euler-Lagrange equations 
we derive the term 


„ p''{x) 2 s'(x)s"(x)p'(x) 


Ho{ sHx)- 


i{x) 


+ 


i{x) 


^ = (^M + m^(x)) g - J^p”ix) + 


(24) 





























from (23) and using Then, from ( |21[ ) we derive the additional tension in the cable h{p) produced 
by its displacement p: 


Hp) = -j^ 


rL J 


s'{x)p'{x] 

i{x) 


dx 


s"{x) 

i{xY 


(25) 


Remark 2. As already mentioned, in classical engineering literature mm one usually replaces 
= -x/l + 'S'( with 1 and neglects the mass of the cable (m = 0) so that one finds (IM, that is, 
s"{x) = with these rough approximations, one can easily obtain also the second order expansion 
of h{p), that is, 

p{x)‘^ dx^ s"(x) + s'(x)p'(x) dx"j (p"(x) — s"(x) — 3s'(x)s"(x)p'(x)) 

= f [ p{x)‘^dx\ s\x) — f f p(^x)dx\ {p"{x) — s"{x) — 3s'{x)s"{x)p'{x)) 

2Tc \Jq J Lc Hq \Jq J 

where we only see the terms of the expansion up to order 2. However, in literature one finds several 
further simplifications. Timoshenko [2T] (see also [221 Chapter 11]) manipulates the functional r(p) 
and reaches the expression 




L, 


\Ho 


p{x) + ^ ) dx ) {s"{x) —p"{x)) 


see [221 (11.16)]. This expression fails to consider some second order terms but inserts some third 
order terms. Further simplifications can be found in literature; Biot-von Karman [71 (5.14)] neglect 
the second term and simply obtain 


AEMg 

t:ie 


y p{x) dx^ {s”{x) — p"{x)). 


This clearly simplifies many computations but this approximation seems not to have a sound justifi¬ 
cation. This is one further reason to stick to a first order approximation of h{p). □ 


Let F be the force in (14), then by taking into account all the energy variations previously found 


we obtain the following equations of motion: 
m^{x)p = 


Hq ... . 2Hr) s'(x)s" ix) ,, , ,, , 

P {x) -- -p (x) - h{p) F{w -p)-Mg, 


^{xfi ^(x 

Miv = —FIw"" — F{w—p) + Mg. 


(26) 

(27) 


We emphasize that the term mg^{x) has disappeared in (26) due to its appearance with the opposite 
sign in (24) and in the gravitational energy. 


2.3 A model for the suspension bridge 

In this subsection we consider the full suspension bridge. We view the deck as a degenerate plate, 
namely a beam representing the midline of the roadway with cross sections which are free to rotate 
around the beam. The midline corresponds to the barycenter of the cross sections, which are seen 
as rods that can rotate around their barycenter, the angle of rotation with respect to the horizontal 
position is denoted by 9. The endpoints of the cross sections are the edges of the plate and they are 
connected to the sustaining cables through the hangers. Then, for small 9, the positions of the two 
free edges of the deck are given by 


w = y ± i sin 9 Ri y ± W . 
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Each free edge of the deck is connected to a sustaining cable through hangers. We derive here the 
equation of this cable-hangers-deck system and, in order to derive the Euler-Lagrange equations, we 
start from the cable-hangers-beam model developed in the previous subsections. 

We denote by pi{x,t) and p 2 {x,t) the displacements of the cables from their equilibrium position 
s(x), see 0: hence, Pi{x, t) — s{x) denotes the actual position of the i-th cable {i = 1, 2). We assume 
that the roadway has length L and width 2i with 2i L. 

• Kinetic energy of the deck. This energy is composed by two terms, the first corresponding to 
the kinetic energy of the barycenter of the cross section, see (10), the second corresponding to the 
kinetic energy of the torsional angle. The kinetic energy of a rotating object is \ where J is the 
moment of inertia and 6 is the angular velocity. The moment of inertia of a rod of length 2i about the 
perpendicular axis through its center is given by ^Mi'^dx where Mdx is the mass of the rod. Hence, 
the kinetic energy of a rod having mass Mdx and half-length i, rotating about its center with angular 
velocity 9, is given by Therefore, the kinetic energy of the deck is given by 


^kd 




dx. 


(28) 


Note that now each cable sustains the weight of half deck, so the boundary value problem (Q for s(x) 
becomes instead 

Hos”{x) = ^ Y + g , s(0) = s{L) = sq . 

• Kinetic energy of the cables. We adopt for both cables and we get 

rL 


m 


Eke = y + P-i) Cix) dx . 


(29) 


• Stiffening energy of the deck. It is composed by two terms, the bending energy of the central 
beam and the torsional energy. As in (|17[), we have 


Eb = ^[ {y"?dx, 


(30) 


whereas the torsional stiffness of the deck is computed in terms of the derivative of the torsional angle: 


E'j' = 


GK 


{e'Ydx. 


2 7o 

Gravitational energy of the deck and cables. These are readily computed and read 


— Mg / ydx and 

Jo 


mg / {Pi +P2)^ix)dx. 

Jo 


Elastic energy in the hangers. Following (15), the energies of the two rows of hangers are 
Ehiiv + JO - Pi) = ^ J K{x)(^(y{x,t) + ie{x,t) - pi{x,t) + s{x) - \{x)'j ^ dx , 

Eh^iv - JO - p2) = ^ J K{x)(^(^y{x,t) - J 0 {x,t) - p2{x,t) + s{x) - X{x)^ ^ dx . 
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Since the weight Mg oi the deck is now supported by two cables, each one supporting Mg/2, the force 
F in (14) acts now on the two edges of the deck and this gives rise to the two forces: 


F{y + ie-pi) 


F{y - £9 - P 2 ) 


k{x) [y{x, t) + W{x, t) - Pi (x, t) + s{x) - A(x)^ 
Kix)(yix,t) + W{x,t) - pi{x,t) + , 

k{x) [y{x, t) - £ 9 {x, t) - P2{x, t) + s{x) - A(x)^ 
Kix)(yix,t) - £e{x,t) - p2ix,t) + , 


with F{0) = Mg/2. Moreover, instead of (13) we have 


2 


k{x) — A ( x )^ 


Ko(^s(x) - X(x)^ 


A(x) 


The Euler-Lagrange equations. Let h{p) be as in (25), then the equations of motion are 

2Hds'{x)s''{x) , 


cl \ " II 

m^{x)pi = -z^—^Pi 

C{xr 

c( ^ " -^0 // 

rai{x)p2 = -rrv?P2 


2Hqs'{x)s"{x) , 


p\-h{pi)+F{y + £9-pi)-^ 


Mg 


P 2 - Hp 2) + F{y - £9- P 2 ) —Y 


i{xY 

Mij = -Ely"" - F{y + £9 - pi) - F{y - £9 - ps) + Mg 
^£^9 = GK9" + £F{y - £9 - P 2 ) - £F{y + £9 - pi). 

O 


(31) 

(32) 

(33) 

(34) 


• Boundary conditions. The degenerate plate is assumed to be hinged between the two towers 
while it is clear that the cross sections between the towers are fixed and cannot rotate. This results 
in the boundary conditions 

y{0,t) = y{L,t) = y"{0,t) = y"{L,t) = 9{0,t) = 9{L,t) = 0 Vt > 0 . (35) 

Since Pi — s denotes the actual position of the cables, we also have 

Pi{0,t) = pi{L,t) = 0 Vt>0, (i = l,2). (36) 


3 Existence and uniqueness for the suspension bridge system 

3.1 Definition of solution 

Up to scaling, we may assume that L = ir: this will simplify the spectral decomposition of the 
differential operator. We endow the Hilbert spaces L^(0, vr), Hq{0, vr), (0, vr) with, respectively, 

the scalar products 

CTT CTT CTT 

{u,v)2 = uv , {u,v)hi = / uv' , {u,v)h2 = / u"v" . 

Jo Jo Jo 

Then an orthogonal basis in these three spaces is given by {ek}'^^i, where 












We denote by H*{0, vr) the dual space of H ^fd(0, tt) and by (•, •)* the duality pairing, whereas we 
denote by the dual space of HQ{0,7r) and by (•, •)! the duality pairing. 

Let ^ be as in Q; we also consider the Sturm-Liouville eigenvalue problem 


Ho , 




= X^{x)u in (0, tt) , u(0) = u{tt) = 0 . 


(37) 


We say that A is an eigenvalue of (37) if there exists u ^ 0 (eigenfunction) such that 


. ° uv' = A f ^{x)uv 
Jo 


It is well-known (see e.g. [U §2.4]) that (37) admits a sequence of positive eigenvalues {Afc} which 
diverges to infinity. Each eigenvalue has multiplicity 1 and the sequence of eigenfunctions {uk} is an 
orthogonal basis of L^(0,7r) and of 77^(0, vr) endowed with the scalar products 

{u,v)^= / i{x)uv Vu, u G L|(0, vr), {u,v)ui = / ° u'v' Vtt, u G 77|(0, vr). 

Jo ^ Jo s(a^) 

Let us emphasize that the corresponding norms || • ||^ and || • are equivalent, respectively, to the 
norms || • ||2 and || • but, for the sake of clarity, we maintain the different notations L| and H^. In 
the sequel, we assume that {uk} is normalized in L|, that is, 

= 1) ll^fcll^i = \/Afc ■ 

We denote by 77^(0, vr) the dual space of 77^(0, vr) and by (•, •)^ the corresponding duality pairing. 

We set 

Mg 


<I>(s) := F(s) — 77^ and 'k(s) = f ^{a)da 

2 Jo 


so that 


$(0) = 4'(0) = 0 


$(s) > 0, ^{s) >0 Vs G 


In the above functional-analytic setting and with this simplification, the equations (31)-(32)-(33) 


(34) may be rewritten as 


m^{x)pi = 

m^{x)p2 = 
Mij 


^(xpPl 

' , AE 
Lc 

■ r-K s'(x)p'Px) 

Jo i{x) 

Hq / 

' , AE 

j-TT S'ix)p'2ix) 

Jo i(x) 


j^ + Hv + JO-pi) 
+ ^{y - le - P2) 




= -Ely”” - $(y + ie- pi) - $(v/ -W- P2) 
= GKO” + m{y -£6- P 2 ) - £^{y + £9- pi) 


(38) 


where ^{x) and s{x) are linked through We add the boundary conditions ([35|)-([^ which we 
rewrite for L = n: 


y(0, t) = v/(vr, t) = y”{0, t) = y"(vr, t) = pi{0, t) = pi^n, t) = (9(0, t) = (9(vr, t) = 0 Vt > 0 . (39) 

We fix some initial data at time t = 0, 


y{x,0) = y^{x), pi{x,0) = p^{x), 9(x,0) = 9°(x} VxG(0,vr) 

y(x,0) = y^(x), pi(x,0) =pI(x), 9(x,0) = 9^(x) VxG(0,vr), 


(40) 
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with the following regularity 


e niJo(0,7r), 0°,^° G iJo(0,7r), 6'^p■ E L2(0, vr). 

We say that {pi,p2,y,0) is a weak solution of (p8|) if 


(41) 


e,pi E ([0, T]; FoHo, tt)) n ([0, T]-L\0, vr)) n ([0, T]-H-\0, vr)) 
y E C0([0, T]; ^2 n ^r)) n ([0, T]-,L\0, tt)) n T]-,H*{0, vr)) 


(42) 


and if the following equations are satisfied 
m{pi,vi)^ + {pi,vi)hi = 4 f [/J" 

IT S'{x)pi2ix) 


m{p2,V2)^ + {P2,V2 )h} = Tf 


/c 


0 €(x) 


2 + {HV + - pi),Vl'j ^ 

'^,V2) +U{y-ie-p2),V2'' 


M{y,w)^ + EI{y,w)H 2 = - (^^{y + W - pi) + ^{y - W - p 2 ),wy 


M 


f{B,V3)i + GK{e,vs)Hi = 


- £9- P 2 ) - £^{y + 16 - pi), U 3 


for all Vi E Hq{0,tt), w E H and t > 0. In the sequel we denote by Xt the functional 

space of solutions: 

{pi,P2,y,9) E Xt <s=> @ holds. 


(43) 


Note that this space already includes the boundary conditions (39): hence, from now on we will not 


mention further (39). 


In view of the energy balance performed in Section 2.3, the conserved (and approximated at second 
er) toti 

m = 


order) total energy of any solution {pi,P2,y,9) E Xt of (38) is given by 

'M „oAo M .0 m, .0 .o, A , f El 

/o\u ^ ^ / J 0 V 2 

s'p[ 


^202 -2 ]dx+ ^(y")^ + 1 dx 


+ 


AE 


2Lr 


m 

6 “ “ ' 2 " ' ^ 

^ AE 

dx + 


s'p '2 


dx \ + 


GK 

ip'lf + {P'2f 


) ^(x) J 2 Lc \Jo ?(x) J Jo \ 2 i{x 

-mg j {pi + p 2 )^{x) dx + j {y + i9 - pi) +'i’{y - £6 - p 2 )^ dx 


2 - s' {p'l + p' 2 ) j dx 

(44) 


/o JO 

where, in the first line we see the total kinetic energy of the bridge and the elastic energy of the 


deck whereas in the second line we see the stretching energy of the two cables. The third line in (44) 


deserves a particular attention: we see there the gravitational energies of the cables and the deck, the 


latter being included in T while the former cancels in the equations (38) due to its presence with an 


opposite sign in both (24) and in the gravitational energy. 


3.2 Existence and uniqueness for a related linear problem 

We consider here the linear and decoupled problem 


m^{x)pi = 
m^{x)p2 = 
My 




+ gi{x,t) 

+ g2{xA) 


fee 


= -Ely'”' + gofyX^t) 
= GKO" + gi{x, t) 


X E (0, vr), t > 0 . 


(45) 
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where gj G C*^([0, vr] x [0,r]) for j = 1, To (45) we associate the initial conditions (40) with the 
regularity as in (41). We say that {pi,P2,y,d) G Xt is a weak solution of (45) if 


Vu G H^{0, tt) , yweH^n H^{0, vr), t > 0 . 


m{Pi,v)^ + {pi,v)jji = {gi,v)2 

m{p2, v)^ + {p2,v)hi = {92,v)2 
M{y,w)^ +EI{y,w)H 2 = {g 3 ,w )2 
ff{e,v)i + GK{e,v)Hi = ig^,v)2 

The purpose of the present section is to prove the following statement. 

Theorem 1. Let T > 0 and let gj G (^^([OjTr] x [0,T]) for j = 1,...,4. For all y^,d^,p^,y^,9^,pj 
satisfying (41) there exists a unique weak solution {pi,P 2 ,y,S) G Xt of (45) satisfying (40). 

The proof of Theorem uses a Galerkin procedure and is divided in several steps; it allows 
to emphasize the delicate role of all the spaces, norms, and scalar products dehned in Section |3.1[ 
Moreover, the proof provides the underlying idea for a numerical approximation of the dynamics with 
a hnite number of modes. 

Step 1. We construct a sequence of solutions of approximated problems in finite dimensional spaces. 
We consider the two basis defined in Section 3.1 and, for any n > 1, we put 

En ■■= spanjei ,... ,en} , Un ■= spanjui,.. .,Un} . 

For any n > 1 we also put 


{Pi)n ■= '^{ptuk)^Uk = - '^\^^^{p'^i,Uk)HlUk 
k=l k=l 

n n 

Vn ■= efc = ^ ^“^( 2 /°, ek)H2 ek , 

k=l k=l 

n n 

■ = efc)2 efc = - ^ ek)H^ ek , 


90 

■'n 


{Pi)n ■= '^{p},Uk)^Uk, 
k=l 

n 

Pn ■= '^iy^,ek)2ek , 

k=l 

n 


k=l 


k=l 


k=l 


SO that 


,,o . „,0 


7/° in , 6*° ^ 6*° , (p°)n ^ Pi in , y^y^, 9^9^, {p])n ^ p\ in 


(46) 


as n —)■ 00 . For any n > 1 we seek {{pi)n, {P 2 )n, Pn, On) such that 

n n 

{Pi)n{x, t) = '^{p.i)\{t)Uk , yn{x, i) = ^ yn(^)efc , 


9n{x,t) = '^9^{t)ek 


k=l 


k=l 


k=l 


and solving the variational problem 


m{{pi)n,v)^ + {{pi)n,v)Hl = {gi,v)2 
m{{p2)n, V)^ + {{P2)n, v)h1 = (^2, v)2 
M{yn,w)2 +EI{yn,w)H2 = {g3,w)2 

^f{9n,w)2 + GK{9n,w)Hl = ig4,w)2 


yv e Un, Vrc ^ En, t > 0 . 


(47) 
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By making n tests on each equation (for v = ui, ...,Un and w = ei, ..., 6 ^) we obtain the 4n linear 
equations 

✓ 

^{Pl)n + ^k{Pl)n = i9l,Uk)2 
m{p2)n + Afc(p2)n = {92,Uk)2 
My^ + Elk'^y^ = ( 53 , 6^)2 

feet + GKke':, = {gi,ek )2 


Mk = 1 ,n . 


(48) 


It is a classical result from the theory of linear ODE’s that this finite-dimensional linear system, 
together with the initial conditions 

(p*)n(0) = {P^,Uk)i, (Pi)n(O) = {p],Uk)^, y^(0) = (y°,efc) 2 , y^(0) = (y\efc) 2 , 

0 ^(O) = (0°,efc)2, ei{^) = {e\ek)2, 

admits a unique global solution defined for all t > 0. 

Step 2. In this step we prove a local uniform bound for the sequence {((pi)n, (P 2 )n, yn, ^n)}- 

We fix some (finite) T > 0; in what follows the Cj’s denote positive constants independent of n, 
possibly depending on T. We take u; = in (47)3 and we integrate first over (0, vr) and then over 
(0, t) for some t G (0, T) to obtain 

M\\yn{t)\\l +EI\\yn{t)\\‘\j 2 = ci + 2 [ {g3,yn)2<ci + 2\\g3\\oo [ ||yn(r)||idr 

Jo Jo 

< C 1 + C 2 f ||yn(r)||idT (49) 

Jo 

where || • ||oo denotes the L°°([0,7r] x [0,T])-norm whereas || • ||i denotes the L^(0,7r)-norm. Here, 
Cl := sup„(M||?/i III -b EI\\y^J\]j2) < 00 . By using the Holder inequality, we see that (^49| implies 


||yn(i)||i < C3 -bC4 l|yn(r)|||dT^ 


1/2 


(50) 


In turn, (50) may be written as 

f{t)<Ci + Ci^/J(^ VtE(0,T), f{t):=[ ||yn(r)|||dr. 


We now recall a Gronwall-type lemma which can be deduced (e.g.) from [3] and [U Lemma A.5/p.157]. 

Lemma 1. Let f G C^(M+) be such that /(O) = 0, 0 < f'{t) < Ci + C 2 \lf{t) for all t > 0 (for some 
Cl, 6*2 > 0/. Then 

fit) + iCi + C2)t Vf>0. 

By applying Lemmato (50) and going back to (49), we obtain constants csjCg > 0 (independent 
of n) such that 

I 2 + l|yn(i) II //2 < C 5 -b cg VtG(0,T). (51) 


Similarly, take v = {pi)n in ([47|)i or n = ip 2 )n in ( 47 ) 2 , to obtain constants 07 , 03 , 09,010 > 0 
(independent of n) such that 


Pi) 


+ Il(yi)n(t)lli7i ^ Cj + Cgt Vf G (0, T) , 
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\\{P2)n{t)\\l + \\{p2)n{t)\\Hi < Cg + Cio ^ Vt E (0, T) . 

In order to obtain these inequalities one also needs to combine the Holder inequality with the equiva¬ 
lence of the norms in and Lg. 


Finally, we take tc = in (47)4. Proceeding as above, we obtain constants cii,ci 2 > 0 (indepen¬ 
dent of n) such that 

ll^n(i) II 2 + 11^1 — *"11 Vt E (0, T) . 

Step 3. We show that {((pi)n) iP 2 )n, Un, On)} admits a strongly convergent subsequence in some norm. 

Let us consider in detail the sequence {yn}, tbe other components being similar. The bound in 
(51) suggests to prove strong convergence in the space 

C'0([0,r]; 772(0, vr))nCH[0,r];L2(0,7r)). 


The equations in (48) show that the components do not depend on n, that is, 

n 

yn{x,t) = . 


k=l 


Take some m > n > 1 and define 


t) ■■= ym{x, t) - yn{x, t) = '^ y^{t)ek 

k=n-\-l 


so that, in particular. 


ym,n(a;, 0) = 0) = l/m “ • 

Consider also the Fourier decomposition of the function g^-. 

00 

9^{x,t) = '^gl{t)ek ; 

k=l 

since 53 E C''^([ 0 , 7 r] x [0, T]) C L 2 (( 0 , 7 r) x (0, T)), we know that 

/ 00 \ 1/2 

( ^ “^0 as n —)• 00 inL 2 ( 0 ,T). 

\k=n-\-l / 


(52) 


Rewrite (47)3 with n replaced by m. Then by taking w = ym,n{t) as a test function, by using the 
orthogonality of the system {cfc}, and by integrating over (0,t) we obtain 

-(l^||ym,n(^) II 2 T F//||ym,^(t) ||gg2 = C*m,n T 2 f (f/S , ym,n)2 — C*ni,n T 2 f { ^ ^ g {x'^Ck, ym,n\ dx 

■^0 -^0 \k=n+l Jo 


^ Cm.n + 2 


^ m \ 1/2 

alitf] \\ym,nA\\2dt 

\k=n-\-l / 


(53) 


where Cm,n = M:\Am, - vWl + EI\\y^^ - yAA- ^7 0 as m,n oo; 

combined with ( 52 ) and ( 53 ), this shows that {yn} is a Cauchy sequence in (//^([O, T]; 77 ^ n//i(o, yr)) n 


C^([0, T]; L2(0, vr)). By completeness of these spaces we conclude that 

3y E C0([0, T]; 772 n 770^0, vr)) n Ci([0, T]; L2 (o, vr)) s.t. 

yn—^y in C'°([0,r]; 772 n 77Q(0,7r)) n C'^([0,r]; 7.2(0, vr)) as n —>•-|-oo , 


(54) 
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thereby completing the proof of the claim. 

For the sequences {On}, {(pi)n}, {(P 2 )n)} we may proceed similarly and prove that 

3pi G CO ([0, T]; H{{0, vr)) n C^ ([0, T]; L2 (o, vr)) s.t. 

{Pi)n-^Pi in CO(^[0,T];i7^(0,7r)^ n C^(^[0,r];L2(0,7r)^ as n+cx), 

39 G CO ([0, tt)) n C^ ([0, T]; L2(0, tt)) s.t. 

Bn^e in CO(^[0,r];/7oi(0,7r)) nC^(^[0,T];L2(0,7r)) as n ^+oo . 

Step 4. We take the limit in ( |47[) and we prove Theorem 

By using (54)-(55)-(56) in (47) we see that {pi,p2,y,0j^is a weak solution of (45) satisfying (40), 
with the additional regularity 

e,piGC^[[0,T]-,H-\0,7r)), yGC^[[0,T]-,H*{0,7T)) 


( 55 ) 


(56) 


following from the equations (45). Therefore, for any T > 0 we have proved the existence of a weak 
solution {pi,P2,y,0) G Xt of (45) over the interval (0,T) and satisfying (40). By arbitrariness of 
T > 0 this proves global existence. Finally, arguing by contradiction and assuming the existence of 
two solutions, we subtract the two linear equations for the two solutions and we obtain an homogeneous 
linear problem; the energy conservation then shows that the (nonnegative) energy is always 0, which 
proves that the two solutions are, in fact, the same. This completes the proof of Theorem 

3.3 The existence and uniqneness result 

The purpose of this final subsection is to prove the following statement. 


Theorem 2. For all y^,9^,p^,y^,9^,pj satisfying (41) there exists a unique (global in time) weak 
solution {pi,p2,y,0) G Xoo of (38) satisfying the initial conditions (40). 

Global in time means here that we can take any T > 0, including T = oo; this explains the notation 
Xoo- For the proof we use a fixed point procedure combined with an energy estimate. The first step 
consists in defining the map which will be shown to have a fixed point; in the sequel we emphasize 
the dependence on time of the nonlocal term. 

Lemma 2. Let T > 0 and assume that y^,0^,Pi,y^,6^,pj satisfy (41). For all {qi,q 2 ,z,a) G Xt 
there exists a unique weak solution {pi,P2,V,0) G Xt of the system 


m^(x)pi = 


mf,{x)p2 = 


Hi 


^ Hi 


0 / 
WPi 


0 / 
^P2 


+ 


+ 


AE 

r 1 

s”{x) 

Lc 

[Jo (ix) \ 

C(x)3 

AE 

■ r s'{x)q'2ix,t) ■ 

s”{x) 

Lc 

[Jo ax) ^ 

C(x)3 


A{x) 

My = —Ely”” — ^{z + (.a — qi) — ^{z — ia — 52 ) 


+ ^{z + ia- qi) 
+ ^{z — (.a — (72) 


M 


^^9 = GKO” + m{z - ea - q 2 ) - i<^{z + ia - qi) 


satisfying the initial conditions (40). 
Proof. We set 

/ ^ \ r 

gi{x,t) = — / 

Uo 


s'{x)q[{x,t) 

^(x) 


dx 


s”{x) 

^(x)3 


+ <k(z F ia- qi ), 


17 







































/ , AE 

g2{x,t) = — 

L^O 


s'{x)q'2{x,t) 

i{x) 


dx 


s''{x) 


+ $(z - la - q 2 ), 


^(x)3 

5 r 3 (x, t) = -$(z + e.a- qi) - ^{z -ia- 92) , 
g' 4 (x, t) = ^^{z — ia — (72) — l^{z E ioL — qi). 
Then gj G C°([0,7r] x [0,r]) for j = 1, ...,4 and we apply Theorem]^ 
Denote by Zt the subspace of Xt such that 


□ 


e,Pi G ([0, T]; FoHo, tt)) n ([0, T]; l2(0, tt)) 
y G CO ([0, T]; n ^1(0, vr)) n C^ ([0, T]; L2(o, vr)) 


(57) 


which is a Banach space when endowed with the norm 

|2 . Il„ Il 2 I II • ||2 


\\{Pl^P2,y,0)\\zT II7'iIIl°°(_H'1) + IIPiIIl“(L 2) + l|P2|lioo(j|^l) + 11^211^00(^2 

+ l|y|lL°°(//2) + 112/11^00(^2) + ||6*||^CX)(^1) + 11^11^00(^2) 


where the norms have the usual time-space meaning. For any yO, 6^,p^, y^, 0^,pj satisfying (41) 

we define the convex subset (complete metric space) 


B := {{pi,P2,y,0) G Zt; (|^ holds}. 

If we take {qi,q 2 , z,a) and {pi,p2,y,0) both satisfying ( [io] ), then Lemmadefines a map 

T:B^B, T{qi,q 2 ,z,a) = {pi,p2,y,9). 

We now prove that for sufficiently small T this map is contractive. 

Lemma 3. Let y^,9^,p^,y^,9^,pl satisfy (41). IfT > 0 is sufficiently small, then the map T defined 


(58) 


in (58) is a contraction from B into B . 


Proof. Consider two different {q\,q\,z^,of) and {q\,q 2 , z^, 0 ^) in B and let {p{,P 2 ^y^~ 
T(y}, q^ 2 } 1 for j = 1) 2. Denote by 

Pi=p\-pI, P2=pI-pI, y = - y^ 0 = 9^ - 9“^ . 


We underline that these notations should not be confused with the initial conditions (40). By sub¬ 
tracting the two systems satisfied by {p{,P 2 ^y ^^^rit {pi,p2,y,9) is a weak solution of 


mf{x) Pi = 


Ho / 

Pi 




AE 


Lr 


C L 4 o 


sfx)Uyix,t)-iqiyix,t)] 

C(x) 


i{x 

+^{z^ -b - g}) - 4>(z^ -b io? - ql) 


dx 


s"(x) 

^(x)3 


(59) 


mf,{x)p2 = 


Ho / 

- 2 P 2 


-b 


AE 


Lr 


C L 4 o 


syx)[{qiy{x,t) - {q^y{x,t)] 
C(x) 


i{x 

+4>(z^ - - ql) - ^{z^ - lo? - ql) 


dx 


s"{x) 

^(x)3 


(60) 


My = -Ely'”' - ^{z^ + la^ - g}) - 4>(z^ - - ql) 

P^iz"^ -b ia^ - ql) + ^{z^ - lo? - ql) 


(61) 
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M 


ee = GKO” + - ql) - <^iz^ + ia^ - g})] 

-l[^{z^ - £a^ - ql) - $( 2 ;^ + la^ - qf)] 


( 62 ) 


satisfying ho mog eneous initial conditions. 

Multiply (59) by pi, (60) by p 2 , (61) by y, (62) by 9. Then integrate with respect to x over ( 0 , 7 r) 
and with respect to t over ( 0 ,t); we obtain 


AE 


m\\pi{t)\\i + \\pimj,, = 2 — 


c Jo UO 


s'{x)[{qiy{x,T) - {qfyix,T)] 

^{x) 


dx 


\ Ljo 


" s"(x) . ■ 

-—^pidx 


t pTT 


+2 


[< 1 >( 2 ;^ +l!.a^ — qi) — ^{z^ — qi )]pidxdT 


0 JO 


AE 


tv ce: J ( 


m\\p 2 {t)\\i+\\P 2 mi,. = 2 — 


c Jo uo 

t rTT 


s'{x)[{qiy{x,T) - {qiy{x,T)] 

i{x) 


dx 


J Ljo 


" s"{x) . ■ 

-—^P2dx 

Ax)-^ 


+2 [^{z^ —ia^ —q2) — ^{z^ —ia —q2)]P2dxdT 

Jo Jo 

M\\yit)\\l+EI\\y{t)\\jj 2 = 2 f f [^{z'^+ia^-qf)-<^>{z^+ia^-ql)]ydxdT 

Jo Jo 

ft /•TT 


Mi 


+2 / / [^{z'^—ia'^— q2) — ^{z^—ia^— q2)]ydxdT 
lo Jo 

rt /*7r 


\\9{t)\\l+GK\\d{t)\\jjr =2£ [^{z^-£a^-ql)-^{z‘^-ea^-ql)]edxdT 

n Jo Jo 


t /'IT 


—2£ / / [^{z'^+£a'^ — qi) — ^{z^+ia‘‘ — qi)]9dxdT. 

Jo Jo 


dr 

(63) 
dr 

(64) 

(65) 

( 66 ) 


Our purpose is to add all the above identities and to obtain estimates. By adding the four left 


hand sides of (63)-(64)-(65)-(66), we find 7 > 0 such that 


M 


m||Pi(i)ll| + l|PiWllli+"^l|P2(t)||| + ||P2(t)||ii+M||y(t)||^ + i?/||y(t)||^. + -^^||0(t)||^ + Gi^||0(t)|| 


m 


> 7(iiPi(t)ii| + WPiimiji + mmi + iiP2(t)ii^i + mmi + + \\emi + wemiii) • (67) 

Next, we seek an upper bound for the sum of the four right hand sides. Two kinds of terms 
appear: the ^-terms and the nonlocal terms. For the term containing <h in (63) we remark that, by 
the Lagrange Theorem, 

|$( 2 ;^ -I- £a^ - grj) - <l>(z^ -|- £a^ - qi)\ < i^ix) (\z^ - z'^j + £\a^ - a^\ + ; 

therefore, by using the Holder inequality several times, 

ft nTT 


0 JO 


< 


[4>(2;^ -|- £a^ — gj) — <1>(2;^ -|- £a^ — qy)]pidxdT 

— ^ J J {\z^ — z^\ + £\a^ — cJ‘ \ + \q\ — |pi|dxfir 

c / {\\J - z'^h + \\a^ - a'^h + hl - ql\\2)\\pi{T)\\^dT . 


( 68 ) 
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Similarly, the terms containing $ in (|64|)-(|65|)-(|66|) may be estimated, respectively, as follows 

P2iT)\\^dT, 



<b-terms in (64) < c / ( \\z^ — z ‘^\\2 + ||a^ — a ^||2 + \\q 2 — 


‘h-terms in (65) < c 


- 2:^112 + ||a^ - a^lb + Iki - qlh + 1^2 “ 


(69) 

r)||2dr, (70) 


<h-terms in (66) < c 


2 ;^ - 2;^||2 + ||a^ - a^'lb + Ik! - q'th + Ik! - 92 II 2 ) |k(r)|| 2 dr. (71) 


,1 


,1 „2 


1 „ 2 | 


10 


Then, we upper estimate the two nonlocal terms in (|63[)-(64): 

rs\x)[{qjy{x,T)-{q^nx,T)] 


/O L4o 


^(x) 


-dx 


J L4o 


s"(a;) . , ' 

- ^Pidx 


^{x) 


dr 


[ [ \iQiyix,T) - {qfyix,T)\d3 

Jo Jo 


< c 


f*7r 


J L4o 


\pi\dx 


dr 


< c 


[\\iQ})\r)-{qfy{T)\\2mr)\\^dT. 

Jo 


(72) 


All together, (68)-(69)-([70l)-([7l|)-([7^ prove that 


the sum of all the r.h.s. of (63)-(64)-(65)-(66) 


< c 


'L (ll^^“^^ll2 + ll«^-“^ll2 + |lk-'?ill2 + l|i?2-'?2ll2) (l|Pl(r)lh + l|P2('r)lh + ||y(T)||2 + ||e(T)||2)c('r 

(ll(k)'(T) - (<?i)'(^)ll2 + II(92)'(t) - (92)'(t)I| 2) (llpd-rllh + l|P2(T)||5)dT 
< c\/t!^||z^— z^|Too(^2)+||q:^— o^|hcxj(^2)+||qi—(} i|Too(2,2)+||92~'?2|Il°°(i, 2)) (l|Pl(''')ll| + l|P2('r)||| + ||y(T)||2 + llk''')ll2)'^’' 

+cv 7 ^||(ij};)' —(qf)'IToo (2,2)+ 11 (92)^ “(92)^11 (^||pi('r)||| + ||p2(T)|||)<ir^ 

< cVrlKk - qj,ql -ql,z^ -z^a^ (IIPl(2-)ll| + l|P2(r)||| + ||p(t)||^ + ||0(T)||i)dT 


1/2 


(73) 


By taking the sum of (63)-(64)-(65)-(66), by taking into account the lower bound (67) and the 
upper bound (73), for all t G(0,rj we obtain 

iipi(i)iii + iipi(i)ii!fi + iiP2(oii| + iiP2(oii!fi + iip(i)iii + iip(t)ii?^2 + iikbiii + ii0(oii!fi 


< c\/T||{g( - q^,q^ - ql,z^ - z^,^! - a^)llz^ 


(IIpi(t)II| + IIP2(t)||| + ||p(t)||| + ||e(r)|||)dr 


1/2 


(74) 


In particular, by dropping the potential part, (74) implies that 

rt 


4>\t)<Ky/4^ with (t){t)= I (Iki(r)|||-h |k 2 (T)|||-h |k(T)||iIk(r)||^)dr 


/o 


and K = cVT\\{q\ - ql,ql - ql,z^ - z‘^,a^ - a^)\\zT ■ 

This differential inequality yields \/(j){t) < Kt/2 for all t G [0, T] which, inserted into the right hand 
side of (74), gives 


ikiWiil + ikiWii^i + \\p2mi + \\p2mi1 + wmwl + ut)fH2 + mmi + pmj,. 
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<cT \\{qi-qi,q 2 -q 2 ,z - z ,a -0)11^^ VtE[0,r]. 

By taking the supremum with respect to t G [0,7"] and by replacing pi, y, 9, this finally yields 

\\{p\ -pl:Pl-phv^ - - 0‘^)\\zt < CT\\{q\ - ql,q\ - ql,z^ - z‘^,a^ - a‘^)\\zT ■ 

Then, for sufficiently small T the map T is contractive. 


□ 


Since a solution of (38) is a fixed point for T (as dehned in (58)), Lemmashows that there exists 
a unique weak solution (pi,P 2 ,y,^) £ Xt of (38) satisfying the initial conditions (40), provided that 
T > 0 is sufficiently small. The conservation of the energy in (44) shows that the solution cannot 
blow up in finite time and therefore the solution is global, that is, T = oo. This completes the proof 
of Theorem O 
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